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Abstract. Nonequilibrium Green's functions represent underutilized means of studying the time evolution of quantum many- 
body systems. In view of a rising computer power, an effort is underway to apply the Green's functions formalism to the 
dynamics of central nuclear reactions. As the first step, mean-field evolution for the density matrix for colliding slabs is 
studied in one dimension. Strategy to extend the dynamics to correlations is described. 
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INTRODUCTION 

Nuclear reactions represent active means of studying nuclei, nuclear interactions and nuclear medium. Central reac- 
tions disturb nuclei more than do the peripheral ones and can lead to nuclear fusion, multifragmentation and dramatic 
compression of nuclear matter, depending on incident energy. Due to the increasing many-body nature of the central 
reactions with higher energies and due to the corresponding population of a multitude of final states, the description of 
those reactions has relied on the time domain and has remained primarily semiclassical. Because of their semiclassical 
nature, the descriptions for central reactions have then remained somewhat disconnected from the descriptions em- 
ployed for nuclear structure, peripheral reactions or giant nuclear excitations. The coarseness of the central-reaction 
approaches, associated to their semiclassical nature, is not easy to quantify and this has proven to be a difficulty in 
linking correctly the observables from the reactions to the properties of nuclear systems under the special conditions 
created in those reactions. 

Historically, the only practical approach to central collisions, that can be considered quantal, has been the Time- 
Dependent Hartree-Fock (TDHF) method ^ and its extensions 10, d]. The starting point for TDHF is the mean-field 
description for nuclei, with individual wave-function components propagated in time according to Schrodinger-like 
single-particle equations of motion. The validity of TDHF relies on the presumption that correlations play a negligible 
role in the dynamics. However, with the weakening of Pauli principle as the excitation energy increases in actual 
reactions, correlations can lead to a fast thermalization of the population of single-particle states as well as to enhanced 
stopping. Such phenomena are completely absent in any TDHF-like description. In view of the issues above, it would 
seem important to develop a quantal approach to central reaction dynamics including the effects of correlations beyond 
the mean-field. These correlations are also known to play a significant role in the initial state |4]. 

The time-dependent Green's functions formalism [5] is a natural candidate for the sought theory. On one hand, equi- 
librium Green's functions can describe the microscopic and macroscopic properties of stationary strongly interacting 
many-body system, accommodating the effects of correlations |4]. On the other hand, an extension of the approach to 
the time-dependent domain, by means of the Kadanoff-Baym (KB) equations, can also incorporate the effect of corre- 
lations onto the time evolution |5]. To date, the application of the nonequilibrium formalism to nonuniform quantum 
many-body systems has been relatively scarce, but this can be in large part attributed to an involved numerical effort 
which should become of a lesser concern as the computational power grows. Besides central reactions, the approach 
could be used to study the response of nuclei to external time-dependent fields [SJ. 

A very important feature of the Green's functions formalism is the fact that, within certain approximations for 
the self-energies, the evolution of the observables is conserving, i.e. the main conservation laws (energy, momentum, 
angular momentum, etc) are preserved during the time evolution |'6','7']. In addition, the propagator formalism naturally 
takes into account the quantal properties of the problem. The non-local features of quantum mechanics are for instance 
reflected into the off-diagonal components of the Green's functions in real space representation. Little is known about 



the structure of these components and their importance in the time-evolution of quantum many-body systems. 



KADANOFF-BAYM EQUATIONS 

The KB equations generally describe the evolution, under rather liberal assumptions, of expectation values of the 
products of two single-particle annihilation and creation operators, at different time arguments, i.e. 2-point Wightman 
functions. In the context of nonequilibrium theory, these are called the single-particle Green's functions. At equal 
time-arguments, the expectation values yield the single-particle density matrix. For realistic modeling of reactions, 
a nonuniform three-dimensional (3D) system needs to be studied. First, however, we shall restrict to the less involved 
one-dimensional (ID) case to get acquainted with the effect that correlations might induce on the time evolution of 
the system. Here, in addition, we shall assume in the following that the initial state of the system is uncorrected, i.e. 
describable by mean-field theory. This approximation will be eventually relaxed, leading to a modification of the KB 
equations which we shall again not treat in here [8]. 

The initial state of the systems is specified in terms of an A-body density operator po at time Iq. The one-body Green's 
functions are then defined as expectation values, with respect to po, of products of Heisenberg-picture creation (x, f ) 
and destruction a{x,t) operators: 

'^^{xi,ti;x2,t2) = i {d'^ {X2:t2)a{xuti)) , (1) 
^>(xi,fi;x2,f2) = ~i{d{xi,ti)a''{x2,t2)) . (2) 

Up to a factor, the time-diagonal Green's function reduces to the one-body density matrix. The KB equations, 
governing the evolution of Green's functions in their arguments, 

dfiE/zf (11)^^(11')+ rdii+(ii)^S(ii')+ r''diE^(ii)^-(ii'), (3) 
dfi^^(ii)E//F(ii')+ rAW+{n)t»{n')+ f''di^^{ii)Y.-{ii'), (4) 

follow from considerations of the equations of motion for creation and destruction operators. In the above, the 
simplified notation 1 = {xi,ti) has been introduced and the retarded and advanced functions are defined according 
to: 

F±(l,2)=F^(l,2)±©[±(ri-r2)][F>(l,2)-F<(l,2)] , (5) 

with standing for a possible singular contribution at f i =t2- A generalized self-energy 1,(1,2) may be introduced fl] 
which accounts for all the interaction effects of the system, including the instantaneous Hartree-Fock contribution, 
E///r(l,2), as well as the ES terms, generated by correlation effects that go beyond the mean-field. 

The complex integro-differential KB equations have to be solved in a self-consistent way, since the self-energies 
depend on the Green's functions that are being solved for. This self-consistency simultaneously guarantees a micro- 
scopic consistency of the theory. In particular, for certain many-body approximations to the self-energy (those of the 
conserving type), one can show that the time evolution induced by the self-consistent KB equations preserves the 
conservation laws obeyed by the system as a whole |l7i[^. These conservation laws can therefore be used as tests for 
the numerical implementation of the equations. 

An attractive feature of the KB equations is their generality. The time evolution induced by Eqs. ^ and (|4|l can easily 
describe different types of correlations, given in terms of different approximations to the self-energy. In addition, the 
equations do not focus on any particular physical system and can therefore be used to study different many-body 
problems. So far, the KB equations have been used to study the time evolution of uniform nuclear matter iSll^l and 
the uniform electron gas [11], as well as inhomogeneous atomic and molecular systems 1 12]. 

The KB equations account for the so-called memory effects. From Eq. (O, one can easily see that that the Green's 
function at times f i and f depends on the Green's functions and self-energies at all the previous times t, f i > f > fo 
and fj/ > f > fo, via the time integrals on the r.h.s of the equation. In consequence, to find a solution of the equations 
one must keep track of all the previous time-steps. A strategy to attack this problem, before addressed in and L13il . 
is discussed further on. 



\ dh 2m dx\\ ^ ' 

^S(ll') = 



_d_ h^d^ 

dtu 2m dx 



v 




-20 -10 10 20 
X [fm] 



-20 -10 10 20 
X [fm] 



-20 -10 10 20 
X [fm] 



FIGURE 1. Time evolution of density for slabs colliding in the mean-field approximation. 



Finally, it is important to note that the KB equations are quantal and, in particular, respect the non-local features 
of quantum mechanics. This is directly reflected in the fact that the evolution, even in the ID case, depends on two 
space and two time variables, (xiJi) and (xii,tii). This doubling of variables might be a disadvantage in terms of 
computational cost, especially compared to other time-dependent approaches which do not share this feature. We 
believe, however, that the increase in computational cost (if any) might be compensated by the advantages of describing 
the system in terms of the Wightman's functions. 



MEAN-FIELD APPROXIMATION 



The rather involved calculations associated with the KB equations simplify when one deals only with the Hartree-Fock 
part of the self-energy. The evolution of the correlation functions is then given in terms of the mean-field and there are 
no memory effects involved in the calculations. The KB equations reduce in that case to the TDHF equations for the 
evolution of single-particle density matrix for the system. This simple framework can be used as a starting point for the 
studies of a time-dependent approach formulated in terms of the density matrix, rather than in terms of single-particle 
wave-functions, as would have been standard for a TDHF approach. 

Since this is meant as a first qualitative example, we have chosen a rather simple local density-dependent mean-field 
of the Skyrme-type: 
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(6) 



with the parameters fo, ti, and a fitted to reproduce the saturation properties of infinite nuclear matter. The local 
density n{x,t) corresponds actually to a 3D system and it is related to the ID density by a scaling factor, n{x,t) — 
—it,'i^{x,t\x\i = xJii = t), where t, = (4no)'''^ "o — 0.16fm^^ is the density of normal nuclear matter. In the 
above, we assume a uniformity of the 3D system in the perpendicular y and z directions, as well as a spin-isospin 
degeneracy of v = 4. 

Since the Hartree-Fock interaction is local, one can implement the time evolution in a rather straightforward way by 
means of the so-called Split Operator Method (SOM) 1 14]. This is based on the fact that, for infinitesimally short time- 
steps Af, the evolution of the density matrix is formally given in terms of the single-particle kinetic K and mean-field 
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FIGURE 2. Intensity plots of the real part of 5f ^ as a function of the x and x' variables. The left, center and right panels represent, 
respectively, the initial state of the system; the state at / = 30 fm/c when slabs interpenetrate and the state at f = 60 fm/c when the 
slabs move away from each other. Note that the white tones correspond to positive values of the function and the darker ones (dark 
blue on-line) to negative values. 



U operators as: 

, , , .(i+0)Al , , .{i+U)Ar' 

^(x,f + Af;x',f' + Af') = e"'— «— ^(x,f;x',f')e' — » — . (7) 

. {K+0)Al .UAl -kM -UAl , 

For short times, one can use the decomposition e s = e^e ~s~e^ +0(At ) to divide the exponential into 
mean-field and kinetic terms. Since the mean-field is diagonal in real space and the kinetic part is diagonal in 
momentum space, we can trivially compute the impact of these two exponentials by applying them separately to 
the real and the momentum space density matrices. By using fast Fourier transforms, we can easily switch from one 
representation to the other and therefore implement the time evolution in a rather efficient way. 

It is worthwhile to comment on the preparation of the initial mean-field density matrix. In order to avoid spurious 
oscillations in various variables for slabs representing the colliding nuclei, it is convenient to start the time-dependent 
calculation from a density matrix that is the ground state for the mean-field. To obtain this ground state, we have 
started from a ID harmonic oscillator density matrix embedded in a harmonic oscillator potential. Then, following 
Eq. (|71), we have very slowly switched off the harmonic oscillator potential, while switching on, at the same time, the 
mean-field. When this process is carried on adiabatically, the ground state of the self-consistent mean-field density- 
matrix is reached at the end of the evolution. After this ground state propagator is obtained, we can repeatedly apply 
the algorithm of Eq. (|7]) to the time evolution of the superposition of boosted ground-state density matrices. This in 
turn gives access to the evolution of one-body observables and total energy. 

In Fig. [T] we show the time evolution for a system where two shells of ID harmonic oscillators have been filled 
for each slab. The slabs have been evolved adiabatically to self-consistency and then boosted by transforming the 
density matrices with the operators e*'^^. The momentum P corresponded to the center-of-mass energy of 25 Me V per 
particle. It is observed that the two initially separated density distributions eventually overlap and interpenetrate each 
other As the system flies apart, a low density region is left behind in which some light fragments seem to be forming. 
This could be naively associated with some form of multifragmentation process. 

One can observe, akeady in this very crude model, that changing collision energy leads to the appearance of 
qualitatively different physical processes. At high energies, the slabs interpenetrate each other and break apart, as 
already commented. At low collision energies, one can observe the fusion of the slabs as well as the subsequent 
evolution of the excited state. These results are in qualitative agreement with the findings of the time evolution of ID 
slabs within the TDHF formalism in terms of wave-functions 1 15]. 

The time-diagonal density matrix depends on two position variables, x and x'. In Fig.|2l we show intensity plots for 
the real part of at three different stages of the collision, specifically at times f = 0, 30 and 60 fm/c. These plots 
give an insight on the structure of the Green's functions away from the jc = x' diagonal. Thus, the initial f = density 
matrices show some spread in direction perpendicular to the diagonal, an effect which can be understood in terms of 
the wave-functions in the structure of the matrix. For one, the support of the matrix is constrained by the support of 
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FIGURE 3. Strategy for solving the Kadanoff-Baym equations in the two time arguments. 



the wave-functions. This yields square outlines where the values of the matrix can be significant. Second, the wave- 
functions signs interplay, producing a band structure along the diagonal. In more dimensions, phases would matter. The 
changing signs and/or phases are generally associated with the changing momentum content in the wave-functions. The 
bigger the momentum spread, the more narrow the spread in \x — x' | . The spread remains similar in the slab collision 
when the slabs interpenetrate at f = 30fm/c. However, the situation qualitatively changes after the slabs collide and 
begin to move away from each. On top of the original structure, the support for the matrix expands, as illustrated 
for t ~ 60fm/c. A larger but fainter square of support emerges, with two corners at the location of slabs on the diagonal 
and two fainter cross-corners persisting into the future. While the matrix values still exhibit sharp maxima along the 
diagonal, the wider region of significant values reflects the fragmentation of original wave-functions. Nucleons from 
the originally populated states can end up in either of the residues and in the space in-between. The portions of the 
wave-functions maintain a phase relationship that would be important for the interference if the portions were ever to 
recombine. 

In the late of the nuclear reactions, though, systems progress through expansions and decays and nearly never 
recombine. Thus, following nuclear interpenetration, we expect to be able to neglect the matrix elements far away 
from the x = x' diagonal without affecting the evolution of the elements close to the axis. In fact, in our preliminary 
tests, we find virtually no effect on the evolution of density on the x = x' axis when ignoring elements farther away 
from the axis than those which are significant for a single slab. Ignoring the elements amounts to the assumption of 
decoherence over farther-away distances. 



CORRELATED APPROXIMATION 

When integrating the KB equations beyond the mean-field approximation, the SOM method alone is not sufficient. 
Even in a ID case, the solution of the KB equations becomes a relatively involved process. At present, we are working 
on the numerical implementation of the equations for nuclear reactions, following closely the implementations already 
existing for uniform matter |j^ ITsd and generalizing them to the non uniform case. In this first implementation, we 
will work within the direct Born approximation for the self-energy. This is the simplest conserving approximation that 
includes effects beyond the mean-field. In the momentum representation, the collisional self-energies are: 

T.^{pa,ta;Ph,th) = J j ^ViPa- Pl)y{Ph- P2)'^^iPuta;P2,tb)'n^iPa- Pl,ta;ph- P2,tb) (8) 

where the function 11^ is a convolution integral of two Green's functions: 

^^iPa,ta;ph,th) ^ J J Y^'^^{Phta;P2,tb)^^{P2-Ph,th;Pl-Pci,ta)- (9) 



In those initial studies, we plan to use a gaussian two-body interaction, V{p) = VQe^~^^ , with parameters chosen 
to reproduce semiquantitative^ the free-space nucleon-nucleon (NN) scattering cross section. This should give the 
first hint on the effect that the inclusion of NN collisions in the formalism has on the time evolution of the different 



properties of the system. Since the Born self-energies are given in terms of convolution integrals, they can be efficiently 
computed using fast-Fourier transform techniques. 

To optimally implement the time evolution given by the KB equations, symmetries should be exploited to minimize 
some of the computational expenses involved in the calculations. A particularly useful symmetry is given by the 
relation: 

[F^(l,2)]* =fS(2,1), (10) 

which is fulfilled by both the Green's functions and self-energies. Therefore, one can store each of the Green's 
functions in one of the half -planes of the (fi,f2) variables and access the other half by using Eq. ( fTOl i. This further 
suggests the implementation of the time evolution following the square structure of Fig. [3] Let us assume that, starting 
from Iq, the Green's functions have been determined for T >t\,t2 > fo and let us denote the rh.s. of (O as !')■ 
To the second order in Af, the Green's function values at T + Af can be then found from 

.E{pi)A/ 1 e — Ti — — 1 

^>(pi,r + Af;/72,f2)-e'^^^>(pi,r;p2,f2) + 2 [^^(/'i,7';/'2,f2)+/>(/5i,r+Af;p2,f2)] — (H) 

with e{p\) = h^p\/2m. An analogous equation holds for A combination of Eqs. (O and (|4|l yields an algorithm 
for evolving the functions along the diagonal in Fig. |3] Advancing of the functions needs to be repeated to achieve 
self-consistency as new function values appear on the r.h.s. as well, e.g. in (fTTl i. 

We expect a general drop in function values ^ for increasing difference \t\ — f2|- Similarly, we expect a drop in 
values for increasing \xi — X2I or \pi — p2\, see Fig.|2] limiting the size of region within the Green's function matrix, 
around the diagonal, that needs to be remembered during evolution. The bookkeeping of values is the particularly 
costly aspect of the numerical implementation of the KB equations. 

In ID already, we expect to gain meaningful insights into how the inclusion of correlations makes the dynamics 
different from the mean-field dynamics or from semiclassical dynamics with nucleon-nucleon collisions. The inclusion 
of correlations may cause turn-on problems early on in the evolution and we intend to test the adiabatic switching on 
of those correlations. Useful for the 2D and 3D calculations should be an experience with treating the far-off elements 
of the Green's functions and even the strategies for discretizing the functions. 
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